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We present a calculation of rj/s for the meson gas (zero baryon number), with the viscosity 
computed within unitarized NLO chiral perturbation theory, and confirm the observation that rj/s 
decreases towards the possible phase transition to a quark-gluon plasma/liquid. The value is some- 
what higher than previously estimated in LO xPT. We also examine the case of atomic Argon gas 
to check the discontinuity of rj/s across a first order phase transition. Our results suggest employing 
this dimensionless number, sometimes called KSS number (in analogy with other ratios in fluid me- 
chanics such as Reynolds number or Prandtl number) to pin down the phase transition and critical 
end point to a cross-over in strongly interacting nuclear matter between the hadron gas and quark 
and gluon plasma/liquid. 



INTRODUCTION 

It has been recently pointed out [l[ that the ratio of the 
shear viscosity to entropy density, rj/s, has an extremum 
at a phase transition, based on empirical information for 
several common fluids, and follow-up calculations by us 
Q and other groups [3] have suggested that 77/s in a 
hadron gas does indeed fall slightly with the temperature 
towards the predicted transition to a quark and gluon 
plasma or liquid phase. Renewed interest in this quantity 
arose after the KSS conjecture [|| about a possible lower 
bound 1 / [Ait) (the existence of a bound had already been 
put forward, on the basis of simple physical arguments, 
in Q) and it is the subject of much current research in 
Heavy Ion Collisions. The precise reach of the bound has 
been under recent discussion, 0, @] and there is much 
interest in finding theoretical or laboratory fluids that 
reach the minimum possible value of rj/s [ol. [Tol] . 

There is good hope that rj/s and even 77 by itself can 
be derived from particle and momentum distributions in 
heavy ion collisions 11 , HI, [HI • 

It has been shown through several examples that, em- 
pirically, rj/s seems to have a discontinuity at a first order 
phase transition, but is continuous and has an extremum 
at a second order phase transition or at a crossover. 

Based on lattice data [3] it is believed that the phase 
transition between a gas of hadrons and a quark-gluon 
phase at zero baryon chemical potential is actually a cross 
over. The result of [l[ however presents a clear disconti- 
nuity. This is of course not a serious claim of that paper, 
but simply an artifact of the very crude approximations 
there employed. We here revisit the issue, improving as 
far as feasible on the hadron-side estimate, and further 
motivating the proposed behavior of rj/s. 



INVERSE AMPLITUDE METHOD IN X PT AND 
HADRON PHASE TRANSITION 

We here improve the very rough calculation of [l[ for 
rj/s on the hadron phase. We have calculated in [TU 
the shear viscosity of a meson gas (that is, the hadron 
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FIG. 1: The viscosity over entropy density of a meson gas in 
chiral perturbation theory unitarized by means of the IAM. 
(z represents the relativistic fugacity e /3 ' M ~ m '). 



gas as a function of the temperature and approximate 
meson chemical potentials, at zero baryon chemical po- 
tential). That work employed the Inverse Amplitude 
Method (IAM) [l]| that gives a good fit to the elas- 
tic phase shifts for meson-meson scattering at low mo- 
mentum, respects unitarity, and is consistent with chiral 
perturbation theory at NLO 



17| . The only explicit de- 



grees of freedom are light pseudoscalar mesons (it, K, 77), 
but elastic meson-meson resonances below 1 GeV appear 
through the phase shifts [18|. 

It is an elementary exercise to divide the calculated 
viscosity from that work by the entropy density of the 
free Bose gas, for N species 
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and plot the result in Fig. [T] Incidently, it can be seen 
in the figure that the holographic bound 2 > 1/47T is not 
violated, which had been claimed in the literature (we re- 
ported this in 0) but independently confirmed by [l9j |. 
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FIG. 2: We improve the hadron-side (low T) estimate of 
[l| that showed the jump in the rj/s ratio in the transition 
from the hadron gas to the quark-gluon plasma, substituting 
the Low- Energy- Theorem of those authors (first order chiral 
perturbation theory) by the Inverse Amplitude Method, that 
agrees with Chiral Perturbation Theory at NLO, and satisfies 
elastic unitarity. We confirm the result of those authors, al- 
though the actual numerical value of rj/s is quite different (as 
should be expected from their calculation reaching tempera- 
tures T ~ 150 MeV but with only the first order interaction). 
One should note that, the calculation being performed at zero 
baryon chemical potential, based on lattice data that sug- 
gest a cross-over between the hadron gas and the quark-gluon 
plasma, and from simple phenomenology this would suggest 
that rj/s should be continuous. 

The reason is that in chiral perturbation theory alone 
the cross section grows unchecked, eventually violating 
the unitarity bound, which induces a very small viscosity. 
Of more interest for our discussion in this work is to ex- 
amine the possible behavior across the phase transition. 
We take the simple estimate for rj/s in the quark-gluon 
plasma from [H, but we use our much improved calcula- 
tion for the low-temperature hadron side (those authors 
employ LO chiral perturbation theory without unitariza- 
tion). The result is plotted in Fig. [21 In addition we plot 
also the phase-shift based phenomenological calculation 
of |4|, that is consistent with ours but somewhat smaller. 

The calculation that [l[ reported shows a discon- 
tinuous jump between the QGP and the hadron gas, 
whereas simple-minded non-relativistic phenomenology 
would make us expect a continuous function with a min- 
imum. Our improved hadron calculation still shows a 
discontinuity, although now the jump at the discontinu- 
ity has opposite sign (our viscosity is larger since the 
meson-meson cross section is smaller due to unitarity, 
instead of being an LO-%PT polynomial) . Since our es- 
timate for rj/s is now approximate only because of our 
use of the first order Chapman-Enskog expansion and the 
quantum Boltzmann equation, both of which are reason- 
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FIG. 3: We plot the dependence of rj/s on the quark mass 
from quark-gluon plasma side by adapting the results from 
[2fl ]. Note that given the non-trivial calculation there, we 
have slightly simplified by taking a constant g in the mass 
correction. The band in the figure corresponds to the inter- 
val g 6 [1,2]. We have also taken all quarks of equal mass 
m s — 120MeV as the maximum possible variation. As can 
be seen, the dependence is small and positive, bringing about 
even better agreement with the hadron side Inverse Ampli- 
tude Method evaluation. 

able approximations, we feel further improvement on the 
hadron side will not restore continuity, and future work 
needs to concentrate on evaluating the viscosity from the 
QGP side. 

To calculate the viscosity in a field theory, a possible 
and popular approach is to employ Kubo's formula in 
terms of field correlators. Another method, based on 
the Wigner function, is to write-down the hierarchy of 
BBGKY equations. 

In either case one can perform a low-density expansion, 
leading to the use of the Boltzmann equation. Employ- 
ing this on the hadron-side. as opposed to the full hier- 
archy of BBGKY equations of kinetic theory, presumes 
the "molecular chaos" hypothesis of Boltzmann, which is 
tantamount to neglecting correlations between sucessive 
collisions. This requires the collisions to be well separated 
over the path of the particle, and induces a systematic 
error in the calculation of order 2 = (m^A), where l/m* 
is the typical reach of the strong interaction, and A the 
mean free path (controlled by the density). To keep this 
number below one requires small densities n(T) < 
If we take as a cross-section estimate 100 mbarn we see 
that the criterion is satisfied up to temperatures of order 
140 MeV (where we stop our plot in Fig. 2). 

We have also estimated the change in rj/s caused by a 
small quark mass, by adapting the results of [23]. Those 
authors provide, within a 2PI formalism, the shear viscos- 
ity of the quark and gluon plasma of one fermion species 
as a function of the fermion mass divided by the tern- 
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perature, for fixed coupling constant. Although we are 
employing, as Csernai et al. do, a coupling that runs 
with the scale (the temperature), the mass correction is 
small, so we can take g — 2 as fixed for a quick eye- 
ball estimate. We normalize the viscosity of [20| to the 
value plotted in Figure [2] at zero fermion mass, and then 
allow the fermion mass to vary. The results are now 
plotted in Figure [3l We plot the extreme case of all 
three light quarks equally massive and with mass equal- 
ing m s — 120MeV. As can be seen, the difference to 
the massless case is irrelevant at current precision and 
does not change the fact that we cannot conclude as of 
yet whether the transition between a hadron gas and a 
quark-gluon plasma/liquid has a discontinuity in rj/s or 
not. The reason that the fermion mass is not so relevant 
in the calculation is twofold. First, even at m s , we have 
T/m s > 1 for any value of T > T c . Since kinetic momen- 
tum transport in a gas is dominated by the fraction of 
molecules (here, partons) with the largest energy allowed 
by the Boltzmann tail of the distribution E > T, and the 
momentum varies as 



2E 2 



P = E(l- 



we see that the small parton mass makes just a correction 
to the momentum of each (efficient) carrier. The second 
reason is that the cross-section, in a regime where pertur- 
bation theory is of any use, is weakly dependent on the 
fermion mass, with a slight dependence brought about 
by the logarithmic running of the quark-gluon vertex. 

Still, given the large uncertainties in our knowledge of 
the quark-gluon medium created in heavy- ion collisions, 
that make difficult to match with the hadron side, we also 
study a simple non-relativistic system where the jump in 
rj/s at the phase transition is very clear. 



LIQUID-GAS PHASE TRANSITION IN ATOMIC 
ARGON 



In prior works it has been pointed out that experimen- 
tal data suggest that first order phase transitions present 
a discontinuity in rj/s and second order phase transitions 
(and maybe crossovers) present a minimum. We will ex- 
amine one case a little closer, for a liquid-gas phase tran- 
sition in the atomic Argon gas, where we will calculate 
the ry/s ratio theoretically and compare to data. The 
empirical data that has been brought forward was based 
on atomic Helium and molecular Nitrogen and Water. 
Quantum effects are very strong in the first at low tem- 
peratures where the phase transition occurs, and the later 
have relatively strong interactions. 

Instead we choose Argon due to its sphericity and 
closed-shell atomic structure, that make it a case very 
close to a hard-sphere system. Thus, Argon is the per- 
fect theoretical laboratory, and sufficient data has been 



tabulated due to its use as a cryogenic fluid. 

The gas phase is therefore well described in terms of 
hard-sphere interactions. In elementary kinetic theory 
one neglects any correlation between sucessive scatter- 
ings. The viscosity follows then the formula 



' mT 

I6d 2 V 7T 



(2) 



where d — 3.42 x I0 5 fm is the viscosity diameter of the 
Argon atom and m = 37.3 GeV its mass 241 ] . 

Experimental data is quoted as function of the tem- 
perature for fixed pressure. The particle density is then 
fixed by the equation of state; therefore a chemical po- 
tential needs to be introduced. In order to calculate the 
entropy density we again use Eq. (pj with N = 1 . 

As said, we keep the pressure P constant, and the 
chemical potential \i varies then within the temperature 
range. In order to obtain \i we simply invert (numeri- 
cally) the function P(T — l/(3,m,n) at fixed tempera- 
ture. The expression for the pressure consistent with the 
entropy above is 



P(T, m, /x) = - 



T 
2^2 



dpp 2 log [I 



(3) 



(we have neglected in both cases the effect of the Bosc- 
Einstein condensate since the gas liquefies before this is 
relevant). The problem has then been reduced to com- 
puting the viscosity at the given temperature and chem- 
ical potential, which we do employing our computer pro- 
gram for the meson gas in the Chapman-Enskog approx- 
imation, with minimum modifications. 

We change variables to absorb the scale and make the 
integrand of order I to: 
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Thus, the final expressions for the entropy density and 
pressure from Eqs. (QJ and ([3]) (once integrated by parts) 
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To treat the liquid Ar phase there is not a very rigorous 
theory. This is because in liquids the momentum transfer 
mechanism is quite complex and does involve the inter- 
action between molecules. Here, our choice of a noble 
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gas is of help since long-range interactions are absent. It 
is common to resort to semiempirical formulas with un- 
clear theory support, or work with formal expressions of 
difficult applicability. We compromise by combining the 
Van der Waals equation of state (that ultimately encodes 
the Lennard- Jones theory for the interatomic potential), 
and use the Eyring liquid theory [22| . 

The Eyring theory is a vacancy theory of liquids. Each 
molecule composing the liquid has gas-like degrees of 
freedom when it jumps into a vacant hole, and solid- 
like degrees of freedom when fully surrounded by other 
molecules. 

This model approach yields a partition function Z for 
a one-species liquid (in natural units) 



E S /N A T f V-V s 



(1 - e- s / T Y 



V s 
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(2n) 3 N A 



, (7) 



from which one can derive complete statistical informa- 
tion about the system [25] ■ One can recognize in the 
second line the partition function of a non-relativistic gas 
for the fraction of atoms with gas-like behavior. The first 
line corresponds to the solid- like behavior. The first fac- 
tor is the partition function of a three-dimensional har- 
monic oscillator. The second term is a correction due to 
the translation degree of freedom, by which an atom can 
displace to a neighboring vacancy. 

The shear viscosity, (like Z itself), turns out to be a 
weighted average between the viscosity of solid-like (first 
line) and gas-like (second line) degrees of freedom of the 
liquid's particles: 
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and to complete the model, 9,n,a, a' , k, E s and V s are 
given in Table U for gaseous Argon. Na is Avogadro's 
number. The V s /V solid- like volume fraction controls 
the weighted average. Note that if this ratio approaches 
1, the viscosity diverges as appropriate for a rigid solid. 

m 

The entropy is calculated as usual taking a derivative 
of the Hclmholtz free energy (A = — TlogZ), 
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For our purposes we also need the liquid density which is 
easily estimated by means of the Van der Waals equation 
of state, that is of some applicability in the liquid phase. 



TABLE I: Liquid Argon parameters which a ppe ars in Eqs. 

and @. All these constants are given in [22|. However, 
k has been modified because we use Eq. ((2| instead of the 
formula that appears in [2^] for the hard-sphere gas case. 



Parameter 



Value 



n 

a — a' 
K 

E s 

Vs 



5.17 meV 
10.80 
0.00534 
0.667 
0.082 eV/particle 
4.16 x 10 16 /m 3 /particle 



This equation takes into account the volume excluded by 
the particles and the attraction between them. In the 
simplest form the Van der Waals equation is: 



(n gas + n 2 Uq ^j (1 - n Uq b) = n Uq , 



(10) 



where n gas and nu q are the particle density of gas and 
liquid Argon, respectively; T is the temperature, 2b — 
47rd 3 /3 is the covolume, that is, the excluded volume 
by the particle (here we take d as a mean value of the 
viscosity radius and the gas radius) and a — 27T C /64P C 
is a measure of the particles attraction (T c — 150.87 K, 
P c = 4.898 MPa). Eq. (TIT)]) is a cubic equation in nu q 
which gives reasonably good results despite its simplicity. 
For this reason, we think that it is not necessary to derive 
a new state equation from the Hclmholtz free energy. 

Putting all together we are able to calculate the rj/s 
ratio in both liquid and gas states. The final result is 
plotted in Fig. [3]where a good agreement with the exper- 
imental data of 23(] is shown. One can see how the KSS 
bound is maintained. Moreover, one can observe that for 
the liquid-gas phase transition rj/s presents a minimum 
and discontinuity at the phase transition (below the crit- 
ical pressure, P c ). Above this pressure, a minimum is 
still seen but the function is continuous. 
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FIG. 4: rj/s (a pure number in natural units) for atomic Ar- 
gon in the liquid and gas phases near the phase transition. 
Solid lines correspond to theoretical calculation described in 
the text, dashed lines are the experimental values given from 
[23| . Note that rj/s is quite independent of the pressure in 
the liquid phase, and that the theoretical curves calculated 
from the liquid side and gas side do get closer together with 
increasing pressure, suggesting as the data that indeed, rj/s 
will be continuous in the cross-over regime. 

CONCLUSIONS AND OUTLOOK 

In this article we have argued, in agreement with pre- 
vious authors, how it is likely that rj/s is a reasonable 
derived observable in relativistic heavy ion collisions to 
pin down the phase transition and possible critical end 
point between a hadron gas and the quark and gluon 
plasma/liquid. We have contributed an evaluation of the 
hadron-side rj/s that simultaneously encodes basic theo- 
retical principles such as chiral symmetry and unitarity, 
and simultaneously produces a practical and good fit of 
the pion scattering phase shifts, by means of the Inverse 
Amplitude Method. In so doing we have updated our 
past meson gas work. Our conclusions are in qualitative 
agreement with those of plj . 

Since our lack of understanding of the non- 
perturbative dynamics on the high-T side of the phase 
transition to the quark-gluon phase prevents us from 
matching asymptotic behavior of rj/s at high T with the 
hadron gas, we have studied this KSS number in a related 
Sigma Model. We find numerically, and confirm with an 
analytical estimate, that keeping the s-channel ampli- 
tude one can isolate a minimum, and within reasonable 
calculational uncertainties, this coincides with the known 
phase transition of the model. A complete analysis is to 
be reported elsewhere. 

Since we are not in possession of a good program that 
can proceed to finite baryon density, we leave this for fur- 
ther investigation. Meanwhile we have investigated the 



past observation that in going from a cross-over to a first 
order phase transition, rj/s changes behavior, from hav- 
ing a continuous minimum to presenting a discontinuity. 
We choose, as very apt for theoretical study, atomic Ar- 
gon. Wc employ standard gas kinetic theory above the 
critical temperature and the Eyring theory of liquids in 
the liquid phase. Whereas the discontinuity in rj/s is 
very clear for low pressures, theory and data are close to 
matching (showing continuity) at high pressures where 
a crossover between the two phases is seen in the phase 
diagram. 

The conclusion is that indeed the minimum of the rj/s 
and the temperature of the phase transition might well 
be proportional. Whether the proportionality constant is 
exactly one could only be established by an exact calcu- 
lation of the viscosity which is not theoretically at hand. 

As a consequence, we provide further theory hints to 
the currently proposed method to search for the criti- 
cal end point in hot hadron matter. If, as lattice gauge 
theory suggests, a smooth crossover occurs between the 
hadron phase and the quark-gluon phase, at least un- 
der the conditions in the Relativistic Heavy Ion Collider 
where the baryon number is small at small rapidity, then 
one expects to see a minimum of viscosity over entropy 
density. In the FAIR experimental program however it 
might be possible to reach the critical end point given 
the higher baryon density (since the energy per nucleon 
will be smaller), and whether the phase transition is then 
first or second order can be inferred from the possibility 
of a discontinuity of rj/s. 
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with Jochen Wambach, Juan Maldacena, Dam Son, 
and Tom Cohen. This work has been supported by 
grants FPA 2004-02602, 2005-02327, BSCH-PR34/07- 
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